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Abstract. We use the statistical hadronization model (SHM) to describe hadron multiplicity yields 
and fluctuations. We consider 200 GeV Au-Au collisions, and show that both event averaged yields 



of stable particles and resonances, and event-by-event fluctuation of the K/n ratio can be described 
within the SHM using the same set of thermal parameters, provided that the phase space occupancy 
parameter value is significantly above chemical equilibrium, and the freeze-out temperature is ~ 140 
MeV. We present predictions that allow to test the consistency of our results. 

> 

The statistical hadronization model (SHM) JH 0. HI] has been extensively applied 
to the study of soft particle production in hadronic systems. When it includes the 
full spectrum of hadronic resonances 0], the SHM, with judiciously chosen (fitted) 
parameters, quantitatively describes the abundances of all hadrons produced in heavy 
ion collisions at all considered reaction energies [5]. 

The ability of the SHM to describe not just averages, but event-by-event multiplicity 
fluctuations has not been widely investigated. Event-by-event fluctuations have attracted 
theoretical ||6l|7|,|8l 9] and experimental interest, both as a consistency check for existing 
models laLZD and as a way to search for new physics W- 

The objective of this work is to determine weather the SHM can describe both yields 
and fluctuations with the same parameters. We obtain a good fit to 200 GeV RHIC 
experimental data including both yields and fluctuations measurements, discuss the 
results in the context of the bulk properties of the matter created at RHIC, and present 
predictions allowing further tests of the model. 

The statistical hadronization model assumes that final states are produced in propor- 
tion to their phase space size. The first and second cumulants of this probability distri- 
bution give, respectively, the average value over all events its event-by-event fluctuation. 

In this work we use the Grand Canonical (GC) ensemble, as implemented in open- 
source software HlQfl . to calculate fluctuations and yields. We motivate this choice by the 
fact that the considered RHIC experiments observe the mid-rapidity slice of the system, 
comprising roughly 1/8 of the total hadron multiplicity. The further assumption of 
(approximate) boost invariance at mid-rapidity allows to image this rapidity slice into a 
domain in configuration space. The matter content in this space domain is expected to be 
in contact and exchanging energy and conserved quantum numbers with the unobserved 
regions. This than creates the GC system we consider on an event-by-event basis. 
If the freeze-out temperature throughout the system is the same, a rather simple model 



can be used to obtain both yields and fluctuations | TH 12]. However, one should note 



that even if practically all produced particles originate in statistical model processes, 
their fluctuations could comprise novel creation mechanisms, related to their formation 
dynamics. However, such novel mechanisms are not present in all observables, and the 
observables we consider here seem to follow non-equilibrium SHM calculations. 

The final state particle yield can then be computed as a function of the particle prop- 
erties, resonance decay tree, freeze-out temperature and fugacities (technical details are 
found in our recent report [13]). While the temperature controls the particle yield depen- 
dence on the mass, the fugacity A describes both the yield of conserved quantities (such 
as baryon number, charge and strangeness) across all particles, and the absolute yields 
which depend on the degree of chemical equilibration. It is common to introduce chem- 
ical potentials associated with each conservation law, /i = Tin A, while the fugacities 
associated with the chemical nonequilibrium condition are called y. 

Detailed balance requires that the particle fugacity be given by the conserved charge 
fugacity to the power of the particle's 'charge', generalized to contain all conserved 
quantities (electrical charge, baryon number etc.) . Thus, for a particle with q, (q) light 
(anti)quarks, s, (s) strange (anti)quarks and isospin I3 we have 

X* = Xf (1) 

However, the condition of chemical equilibrium might no longer hold when the fireball 
is rapidly expanding. Thus, chemical parameters acquire a kinetic (time-dependent) 
component parametrized in terms of phase space occupancies, 

h = K q f q + H + \ (2) 

where y q = 1 , y s = 1 at equilibrium. Even in chemical nonequilibrium the particle fugac- 
ity Xi is the parameter controlling the particle yield, and the first and second cumulants 
can be calculated from the partition function in the usual way ifLUflsll . 

If the expanding system undergoes a fast phase transition from a Quark Gluon Plasma 
(QGP) to a hadron gas (HG), chemical non-equilibrium ifl^l and super-cooling lfl7ll 
are expected to arise given the requirement that entropy has to increase while the 
transition occurs from a high to low entropy density phase. The virtue of a hadronization 
temperature near 140 MeV, and an over- saturated phase space (y q ~ 1.5, y s ~ 2), is a 
match of both the energy and entropy density between the QGP and HG phases. 

Fits to experimental data at both SPS and RHIC energies support these values of 
y q)S . Moreover, best fit y qjS > 1 arises for a critical reaction energy ^ (corresponding 
to the energy of the K/n "horn" EH) and system size II 1911 . as expected from the 
interpretation of y„ as a manifestation of a phase transition. However, even though the 
fits performed in [5] strongly favor the chemical non-equilibrium, they do not rule out 
equilibrium models. The equilibrium model remains marginally compatible with data 
giving a less convincing but statistically still non-absurd validity.This happens since 
in a fit considering only particle yields, the chemical non-equilibrium phase space 
occupancies y s and y q correlate with freeze-out temperature [5], making a distinction 
between a higher temperature equilibrated freeze-out (y q = l,y s < 1) scenario and a 
supercooled scenario where y q>s > 1 difficult. When full chemical nonequilibrium is 
allowed for, y q — 1.6 and T = 140 MeV is found. The best fit freeze-out temperature 



when full chemical equilibrium is assumed varies between studies, ranging from T = 
155 MeV in latest SHARE based studies [5], T = 165 MeV for those carried out 2 years 
ago by STAR experimental group [ 20] , to T = 111 MeV offered in the initial RHIC data 
exploration in which strange particles were not yet fully incorporated 12111 . 

The study presented in II 1 311 makes it clear that the dependence of the fluctuation 



on T and y q is different, allowing us to independently determine these two variables. A 
higher temperature decreases the fluctuations with respect to the Poisson value a = 1, 
expected for a Boltzmann distribution, since it introduces greater particle correlations 
arising from increased resonance decay contribution. Increasing y q rapidly increases 
fluctuations of quantities related to pions, due to the fact that at y q > 1 X K rapidly 
approaches e m *i T , giving fluctuations an extra increase compared to yields mm. 

By virtue of the implied physical picture, equilibrium models generally assume a 
long time span between chemical (particle production) and thermal (particle scattering) 
freeze-out, which would alter considerably the multiplicity of directly detectable reso- 
nances. In the chemical non-equilibrium supercooled freeze-out picture, however, it is 
natural to assume that particle scattering after emission is negligible [ 16] and thus one 
can in most cases assume that the thermal and chemical freeze-out temperatures are the 
same. Hence, a reliable way to probe the re-interaction period would be instrumental for 
our understanding of how the fireball produced in heavy ion collisions breaks up. 

We have recently shown HOI that a comparison of fluctuations to directly detected 
resonances probes the interval between chemical and thermal freeze-out. Consider, for 
example, o K + 1 % - . The numerator and the denominator terms in this ratio are linked by a 

large correlation term due to the A^*°(892) — > K + n decay. This correlation probes the 
^*°(892) abundance at the initial chemical freeze-out, since subsequent rescattering of 
K* — > Kit decay products or on-shell K* regeneration from in-medium K + Tt pairs does 
not alter the final abundance of TZ + and a K . On the other hand, a direct measurement 
of the A^*°(892) /K~ ratio through invariant mass reconstruction measures the ^*°(892) 
abundance at thermal freeze-out, after all rescattering/regeneration ceased. Hence, com- 
paring the K + /n~ fluctuation to the ^*°(892) /K~ ratio provides a gauge for the effect 
of the hadronic reinteraction period on particle abundances. A strong constraint arises in 
a model where chemical and thermal freeze-out coincide, as both observables must be 
described with the same set of statistical parameters determined by global yields of other 
particles. In this way one can argue that the K + /n~ fluctuation and the K*°(892)/K~ 
relative yield measured by invariant mass method offer a decisive test of the chemical 
non-equilibrium sudden hadronization reaction picture. 

While they are phenomenologically powerful, fluctuation measurements are also vul- 
nerable to systematic effects which need to be carefully considered. Volume and cen- 
trality fluctuations, difficult to describe in a model-independent way, can be taken care 
of by considering event by event fluctuations of particle ratios, where the fluctuation in 
volume cancels out, event by event, to first order. This leaves, however, possibly large 
effects due to limited experimental acceptance. These can usually be subtracted by con- 
sidering fluctuations measured within fake events, created using mixed event techniques 
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ll23ll . Such "static" fluctuation, with, by definition, no correlations between particles, can 
be described by a purely Poisson term, where fluctuation is governed by particle yields: 



° stat (Ni) + (N 2 ) (4) 

It can also be seen that mixed-events, made from tracks measured in a given detector, 
also contain the effect due to that detector's acceptance. 

Subtracting o stat from the total fluctuation leaves the "dynamical" fluctuation term: 



d*» = yl (5) 

which comprises the physically interesting effects such as resonance decay correlations, 
Bose-Einstein correlations, and eventual new dynamics. Provided certain assumptions 
for the detector response function hold (see appendix A of I23I1 '). this dynamical fluctu- 
ation is a "robust", detector-independent observable. 

When using measured o dyn s in fits to experimental data, the data-sample should 
include o dyn , particle yields (which are needed to determine the Poisson contribution to 
Ostat as per Eq. @) and particle ratios. We have performed a fit incorporating all STAR 
ratios given in Ref. 112011 . with the exception of the A ++ / p, which the STAR collaboration 
has since begun to reevaluate, and we also for the present ignore the Q./Q. > 1, which 
cannot be fitted with the SHM 1 13, 0]. On the other hand, we have included in the 
fit procedure the preliminary value of o$? n measured by STAR 1E2L as well as the 

published yield for (j) lEll and %~ ll26ll . We assume ll26ll27ll full feed-down correction 
for Ks.l — ► and A — ► % weak decays, and no correction for A — ► p. 

The fit parameters include the overall normalization, the freeze-out temperature T, 
X q>s j 3 and jq.s- We also require, by implementing them as "data-points", strangeness, 
electrical charge and baryon number conservation: (s — s) = (Q) / (B) — 
((Q) j (B))au = 0.4. The fit's statistical significance (Ptrue) profile is shown in the 
left panel of Fig. [TJ To obtain a profile for P true a fit was performed for each fixed 
parameter point on the abscissa, all fit parameters except the one shown on the abscissa 
were varied. It can be seen that the fit tightly constrains y q well above the equilibrium 
value, accompanied by T ~ 140 MeV, in good agreement of the prediction of the 
supercooled hadronization scenario. 

The right panel of Fig. [l] shows the sensitivity of to y q and temperature, and 
explains why the correlation between T and y q disappears when fluctuations are taken 
into account. As can be seen, a fit assuming the chemical equilibrium (y q = 1, T = 155 

MeV |0]). misses o^Jk by many standard deviations. On the contrary the chemical 
nonequilibrium fit seems to be right where this fluctuation is measured. Introducing 
exact conservation for strangeness within the observed window (canonical ensemble) 
would decrease the theoretical gk 12a], thereby increasing the chemical equilibrium 
theory to experiment discrepancy. It is only through y q > 1 that o d ^J K increases to the 
point where it becomes compatible with the experimental value. 

An eye assessment of the fit's goodness is provided by the left panel of Fig. El As 
can be seen, the fit gives an adequate description of all particle yields, including the 
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Y q T[MeV] Tq T[MeV] 

FIGURE 1. Left: statistical significance (Ptme) profiles for freeze-out Temperature T and light quark 
phase space occupancy y q obtained in the fit shown in Fig. [2] Right: Sensitivity of £7^™. to freeze-out 
temperature and y q for the fit in Fig. [2] At each point in the abscissa, a fit is performed , with only the 
particle yield data-points used, varying all fit parameters except the one on the abscissa. o c ^J % is evaluated 
using the best fit parameters. The dot-dashed lines refer to the experimental limits [22] for the a K + j n + 
(blue, lower values) and O k - / K - (red, higher values). See ItHi for further details. 



resonance (iT (892) + K*°(S92))/K- and A(1520) /A. It can also adequately describe 
the event-by-event fluctuations of K + / n + or K ' j%~ . 

K 1%^ fluctuations do not directly test the simultaneous freeze-out hypothesis, since 
K n and K + iz + are not correlated by resonances. To test sudden freeze-out, we 
have used the best fit parameters to predict the yields of several resonances subject to 
current experimental investigation (p°,/ (950),A ++ ,E* + (1385),E*(1530)) as well as the 
dynamical fluctuation of the ratio of their decay products (n + / ll~ ,pj ' ,A/ ^ ± ,S/^ ± ). 
The result is shown in Fig El right panel. 

Note the significant difference between ratios such as p/n + and p/n~ , while the fluc- 
tuations of p/K + and p/n~ are substantially identical. This systematics, which repeats 
itself when the A/n ratio is considered, is due to the correlations provided by the leading 
resonance decay (A — > pn, E*(1385) — > An). Thus, the combined measurement of the 
resonance and the ratio of decay products yields a very powerful constraint on the simul- 
taneous freeze-out model considered here, and it will be interesting to see to what extent 

will the model agree with data. In particular, the difference between <J^+" n+ and a^, % ^ 
is intriguing, since the isospin chemical potential required to reproduce it (A/3 ~ 0.96) 
is excluded through ratios such as the tz + /iz~ . It remains to be seen weather this result 
is due to experimental systematics not analyzed within the preliminary measurement, or 
weather additional theoretical insights are needed to describe it. 

In conclusion, we have shown that the SHM can describe both the yields and event- 
by-event fluctuations measured so far in RHIC 200 GeV Au-Au collisions, provided that 
phase space is saturated above equilibrium and the system is super-cooled with respect 
to the phase transition temperature. We have justified this scenario in the context of a fast 
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FIGURE 2. Left: Fit of preliminary 200 GeV data, including the ^ ± /^ ± fluctuations and the K* (892) 
and A*(1520) resonance Right: Predictions of resonances and event-by-event fluctuations of their decay 
product ratios with the best fit parameters. 



phase transition from a high-entropy phase, and argued that the simultaneous description 
of yields and fluctuations is consistent with an explosive freeze-out, where interactions 
after hadronization are negligible. We have predicted experimental observables suitable 
for testing this model further, and eagerly await more published data to determine to 
what extent can the SHM account for both yields and fluctuations in light and strange 
hadrons produced in heavy ion collisions. 
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